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Quasi-localization and quasi-mobility edge for light atoms mixed with heavy ones 
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A mixture of light and heavy atoms is considered. We study the kinetics of the light atoms, 
scattered by the heavy ones, the latter undergoing slow diffusive motion. In three-dimensional 
space we claim the existence of a crossover region (in energy), which separates the states of the light 
atoms with fast diffusion and the states with slow diffusion; the latter is determined by the dephasing 
time. For the two dimensional case we have a transition between weak localization, observed when 
the dephasing length is less than the localization length (calculated for static scatterers), and strong 
localization observed in the opposite case. 
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Mixtures of different species of cold atoms present an 
interesting field of many particle physics. Two or more 
different types of atoms can be mixed, where one type 
of atoms can be relatively light (e.g. 6 Li), and the other 
type is heavy (e.g. Rb). Quantum tunneling of light 
atoms is a phenomenon, interesting both from an ex- 
perimental and theoretical point of view. The heavy 
atoms serve as slow moving scatterers for the light atoms. 
Lately, it was realized that ultracold atomic gases ap- 
pear very convenient for experimental studies of Ander- 
son localization of the light atoms, both for the case 
of Bose-Einstein condensates, and for fermionic gases 

Kinetics of classical particles in a disordered medium 
can be described by the Boltzmann equation. The most 
drastic manifestation of the difference between the ki- 
netics of classical particles and that of quantum ones is 
Anderson localization. It is well known that for d = 1 
and d = 2, where d is the dimensionality of space, all the 
states are localized, and for d = 3 there exists a mobility 
edge E c , the energy which separates the states with finite 
diffusion coefficient and states with the diffusion coeffi- 
cient being exactly equal to zero. (For reviews, see e.g. 



Refs. [H, lH, 14 1.) All this is true provided the disorder 
is static. A natural question arises: what happens with 
this picture when the scatterers slowly move? 

To answer this question we need some quantitative the- 
ory of localization. As such we will use the self-consistent 
localization theory by Vollhard and Wolfle [l5| . Of cru- 
cial importance in the above mentioned theory are max- 
imally crossed diagrams (the sum of all such diagrams is 
called the Cooperon) for the two-particle Green function. 
The calculations of these diagrams for the case of moving 
scatterers were done in the paper by Golubentsev [l6j |. 

One should notice that we consider the heavy atoms as 
classical objects whose diffusive motion is not affected by 
localization effects. On the other hand, we consider the 
light atoms as quantum objects. Thus the temperature 
of the atom gases should satisfy the inequalities 
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where M and N are the mass and concentration of heavy 
atoms respectively, m and n are the mass and concentra- 
tion of light atoms and T is the temperature. The large 
ratio between the masses of the two types of atoms con- 
sidered is crucial for the applicability of the methods used 
in this work also because following Ref . [lj| , we shall ig- 
nore the change of the energy of the light atoms as a 
result of a scattering by a heavy one. 

In the first part of the present paper we reproduce 
the results by Golubentsev (trivially generalized for the 
arbitrary dimensionality of space) . In the second part we 
use the results for the Cooperon as an input for the self- 
consistent localization theory, which we modify to take 
into account the slow motion of scatterers. In the third 
part we discus the results obtained. 

The quantum particles are scattered by the potential 

V(r,t) = V^2Hr-r a (t)). (2) 

a 

Define the correlator 

K(r -r',t- t') =< V(r, t)V(r', t') > . (3) 

To leading order in the density of scatterers we have for 
the Fourier component of the correlator 

K(q,t) = V 2 ^J exp{( iq (r-r')} 
xdrdr' (r - r a (i)) £«S (r' - r' a ,(0))\ 

a a' I 

= V 2 J2 < ex P - r o (0))}) - nV 2 f(q, t), (4) 

a 

where n is the scatterer density. We consider the case 
when the scatterers undergo slow diffusive motion. In 
the ballistic case 



f(q,t) = cxp 



In the diffusive case 
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/(q,t) = exp 



V 2 Tt„ 



2M 



\t\ , \t\ » r imp , (6) 
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where we have used the fact that 

2 dT 
< v >= — 

and Timp is the mean free time of the scatterers. 



(7) 



For the Cooperon we get [1( 



f t ,)dt' } dt,(8) 



C E (q) = exp (-D(E)qH - ± jf (1 



where -E is the energy of each of the two quantum particle 
lines in the Cooperon diagram, and q is the sum of their 
momenta (see Fig. 1). Also 



nV 2 k 2 /nv d = 3 
nV 2 k/v d = 2 
nV 2 /v d = 1 



(9) 



We'll assume that t <C n mp . The quantity f t is /(k) 
averaged with respect to the iso-energetic surface. We 
obtain 




|*| < T, 
1*1 > T r , 



imp 



where 



( 2k 2 T 



(10) 



(11) 



For d = 3, y 3 (ac) = (1 - e~ x )/:r 16]. For d = 2 



Using the integral 
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f(k(s-s'),t). 



-f dde- A(1 - cose) =e- A I (A), (13) 
t Jo 

where Jo is the modified Bessel function, we obtain 

V2(x) = e-*l 2 h{x/2). (14) 



For d = 1 



(15) 



Eq. dHJ) can be easily understood if we compare dia- 
grams for the Diffuson (the sum of all ladder diagrams) 
and the Cooperon in Fig. 1. The Diffuson does not have 
any mass because of the Ward identity. In the case of 
the Cooperon, the Ward identity is broken, and the dif- 
ference [1 — /(*)] shows how strongly. The interaction line 
which dresses the single particle propagator is given by 
the static correlator, and the interaction line which con- 
nects two different propagators in a ladder is given by the 
dynamic correlator. The time-reversal invariance in the 
system we are considering is broken due to dephasing; 



2 

(a) 



t-l t-2 

(b) 



FIG. 1: Diagrams for the Diffuson (a) and the Cooperon (b). 
Solid line is dressed quantum particle propagator, dashed line 
connecting points r, t and r', t' corresponds to K(r — r', t— t'). 



the diffusion pole of the particle-particle propagator dis- 
appears, although the particle-hole propagator still has 
a diffusion pole, which is guaranteed by particle number 
conservation. (The weak localization effects for the case 
of inelastic electron-phonon scattering were discussed by 
Afonin et al. [lj|.) 

Considering the limiting cases, from Eq. ([5]) we obtain: 



(i) in the case 2k 2 T/M < t/t^ 



C E (d) = / exp [-D(E)(i 2 t - t 2 /r 2 (E)] dt, (16) 
Jo 



where 



2M t 
k?T : 



imp 



1/2 
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(ii) in the case r/if mp < 2k 2 T/M < 1/r 2 



C E (q) = / exp [~D{E)ci 2 t - t 3 /^(E)] dt, (18) 
Jo 



(12) where 



^ k 2 T J 



1/3 



(19) 



Thus we obtain the crucial parameter - the dephasing 
time Tip. 

The results for the dephasing time (up to a numeri- 
cal factors of order one) can be understood using simple 
qualitative arguments. Consider the ballistic regime. If 
a single collision leads to the quantum particle energy 
change SE, the dephasing time could he obtained using 
Eq. EU 



t v SEJ^ 

T 



2tt, 



(20) 



where t v /t is just the number of scatterings during the 



time T v . So in this case 



{8Ef 



(21) 



If we notice that 1/t\ is the averaged quantum parti- 
cle energy change in a single scattering, SE, we imme- 
diately regain Eq. ([11]). Eqs. (gDJ) and (JSTJl also imply 



3 



that if the scattering is quasi elastic (and slow motion of 
scatterers means just that), the energy relaxation time is 
much larger than the dephasing time [191 ] . Hence we have 
the right to ignore the Doppler caused cumulative energy 
shift, which otherwise would have lead to the appearance 
of the Diffuson mass. 

Inserting Eq. (TT5)) into the self-consistent equation, for 
the diffusion coefficient D we obtain 



Dq{E) 
D(E) 



1 



— Y 

4:ir 2 mk ^— ' 
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C B (q) 



(22) 



where Do is the diffusion coefficient calculated in the 
Born approximation 



Dn = -v 2 t: 



(23) 




FIG. 2: The results of numerical solution of Eq. (|26[l {t v /t = 



10) for g{x) = x 2 (solid line), g(x) 
g(x) = x (dots). 



x (dashed line), and 



v is the particles velocity, and the momentum cut-off 
|q| < l/£ is implied, where I = kr/m is the mean free 
path. Thus we obtain 



Do 
D 



1 



dt 



dq q 



d-l 



ll 



irmk jo 

x exp [-Dq 2 t - g(t/T v )] , (24) 



where g(x) is some function which goes to infinity when 
x goes to infinity as some power of x higher than one (in 
the particular case of ballistic regime g{x) = x 3 , and in 
the diffusive regime g(x) = x 2 .) 

Introducing dimcnsionless variables we obtain 



Do 
D 



1 

1 + - 
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7T {uy- 1 

x exp 



dt I dqq a 
'o 



1 D ~2~ f~ I \ 

a L>o 



(25) 



Thus we have obtained an algebraic equation for D / Do, 
which (equation) depends upon two parameters: t,,/t > 
1 and kl, which can be arbitrary. 

Let us start analysis of this equation with the case 
d = 2. Calculating the integral with respect to q we 
obtain 



D 
Do~ 
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nkl 



r°° dt 


Dt ' 




1 - e 


L t 





-g(tT/T v ) 



(26) 



Let us make the assumption (which we'll justify a poste- 
riori) 



Dt v /D q t > 1. 
To calculate the integral 

-s(t) 



r°° dt 




lo t 





A> 1, 



(27) 



(28) 



let us divide the region of integration by choosing some 
x satisfying 1/A <C x <C 1. We obtain 
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ln(Ax) 


— In a; 


= In A. 



(29) 



(In Eq. we have ignored all numerical factors of 

order 1 in the argument of the logarithms.) Hence, Eq. 
(|26p can be presented in the form 



^=l-J-ln(^ 

Do TTkl \D T 



(30) 



Solution of Eq. ([30)1 is particularly simple in two lim- 
iting cases: l v -C C an< i 3* £i where l v = vt v is the 
dephasing length, and £ = le JL 5~ is the localization length 
[14j . In the former case we obtain just weak localization 
corrections 



^ = l-4-ln^ 
Do irkl t 



and in the latter case 



£, 2 

D = — . 



(31) 



(32) 



We see that in both cases the assumption (|2"7) is satisfied. 

Results of a numerical solution of Eq. (|2"6"]) for g(x) = 
x 3 , g(x) = x 2 and g(x) = x are presented on Fig. 2. One 
can see that the curves for D/D-j are practically indistin- 
guishable. Thus the exact form of the function g(x) is 
not important. All the relevant information is contained 
in the dephasing time, determined by the parameter r^. 

Notice, that the quantum diffusion of particles scat- 
tered by the slow moving scatterers turns out to be sim- 
ilar to the case when there are two separate scattering 



4 



mechanisms: strong elastic scattering causing relaxation 
of momentum, and weak inelastic scattering due to say, 
phonons, causing dephasing (except for the definition of 
t„). Strong dependence of the diffusion coefficient for 
d = 2 upon the ratio of the dephasing and the localiza- 
tion length (for the case of two scattering mechanisms) 
was thoroughly discussed in Refs. @, [2^, 21 , 2p. 

As it was noticed by Gogolin and Zimanvi [20j|. there 
is a lower bound of temperature for the validity of Eq. 
(f32l At low enough temperatures variable range hop- 
ping, which is of course not taken into account by the 
self-consistent localization theory is the main diffusion 
mechanism. So Eq. (f24|) is valid, provided 



(33) 



T» AE, 



where AE is the average energy difference between neigh- 
boring localized states. Eq. ([33|) can be presented as 



^ ml 2 



(34) 



On the other hand, inequality l v 3> £ after substitution 
of Eq. gives 



M 



1 



and after substitution of Eq. p^]) gives 



M 1 



-Tthl 



ml 2 



(35) 



(36) 



We again obtain see the importance of the large param- 
eter M/ra. 

In fact, Eq. (f3"2"|) is valid both for d = 1 and d = 3 
(in the latter case, provided we have localization in the 
absence of dephasing). Taking into account the numer- 
ical results obtained for d = 2, for the purpose of semi- 
quantitative analysis we may approximate Eq. (|24[) by 



1 



nmk 



dt 



i/' 



dq q 



d-l 



exp [-DqH - t/r v ] . 



(37) 



Calculating the integral with respect to t we obtain Eq. 
(13711 in the form 



D 

A) 



1 dqq 



nW-iJo q 2 + 1 



For d = 2 we obtain 



1 



£ = 1— bin 



-Dry 
I 2 
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(38) 



(39) 



which in our approximation coincides with Eq. ([30 
For d = 1 we obtain from Eq. (|37|) 



1 



1 y/~D^p _ x 

— tan — . 

7T / t 



(40) 



Again ignoring numerical multipliers of order 1 we obtain 



D = Dn 



(41) 



If we take into account that for d = 1 we have £ ~ I, we 
see that Eq. (|4"Tj) is equivalent to Eq. (f3"2"|) . One must 
admit, however, that for d — 1 the self-consistent local- 
ization theory should be handled with care. In addition 
interaction between quantum particles, not considered in 
the present paper, may strongly influence the localization 
processes [23j. 

For d = 3 from Eq. (|3T[) we obtain 



D 
A) 



= 1 - 



ir{kl) 2 



1 



I 



: tan 



(42) 



One can see, that for d = 3 (similar to the case d = 2) 
Eq. (f3"2"l) ceases to be valid when the localization length £ 
becomes large enough, which happens when the param- 
eter hi approaches the critical value y/S/n from below. 
In fact, in this region Eq. (|42]) can be presented as 



D 



= 2V3^{\c - A) 



I 



yfD, 



: tan 



(43) 



where A = 1/nkl. and A c = l/y&K. After assuming that 
the term 2-\/37r(A c — A) can be ignored with respect to the 
second term in the rhs of Eq. (|4"3"|) , and that Dt v /1 3> 1 
we obtain 



D 



r 2 



(44) 



Now checking the assumptions and taking into account 
that in the critical region 14| £ = l/\X— A c |, we see that 
Eq. (|33|l is valid, provided 

£ > ^ 2/3 ^ /3 . (45) 

The results of numerical solution of Eq. (|42|) are pre- 
sented on Fig. 3. 

Notice that in accordance with Refs. [13, Ei| the de- 
phasing time dependence of the diffusion coefficient can 
be obtained from its frequency dependence by replacing 
u> by iTm. Eq. (|24|hi the absence of dephasing but for 
finite frequency is [14j 



Do 
D(lu) 



I roo l-l/l 

1 + r dt dq q*- 1 

D{uj)q 2 t + iut)] . 



x exp 

The localization length is defined [1 



as 



£ = lim 



D(u) 



(46) 



(47) 



Analyzing the solution qualitatively, we may substitute 
1/t v for — iu> into the definition of the localization length 
(HrJ)) and obtain Eq. 
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FIG. 3: The results of numerical solution of Eq. (|42|) for 
t^/t — 10 (dashed line), t v /t — 100 (dot-dashed line), and 
t v /t = 1000 (solid line). 



CONCLUSIONS 

We considered the influence of slow random motion of 
random scatterers on the localization of quantum par- 
ticles. It turned out that whenever the states of the 
quantum particles were localized, under the assumption, 
that the same scatterers are static, taking the motion 
of the scatterers into account leads to a finite value of 
the diffusion coefficient. In particular, for the three di- 
mensional case, there exists a narrow crossover region in 
energy space, which separates the states with high and 
low diffusion coefficient, the latter being inversely pro- 
portional to the dephasing time. (For the states with 
fast diffusion the dephasing is irrelevant.) Like the po- 
sition of the mobility edge in the case of static scatter- 
ers, the position of this crossover region is determined by 
the criterion that the mean free path is of the order of 
the quantum particle wavelength. This crossover region 
we call the quasi-mobility edge, and the phenomena in 
general we call quasi-localization. For the two dimen- 
sional case we have a transition between weak localiza- 
tion, observed when the dephasing length is less than the 
localization length (calculated for static scatterers), and 
strong localization observed in the opposite case. 

The main application of our results we see as lying in 
the description of kinetics of ultracold gases. However, we 
would like to mention possible application of these results 
to at least one other field. In our previous publication 
[iij ] , we studied the influence of dephasing on the Ander- 
son localization of the electrons in magnetic semiconduc- 
tors, driven by spin fluctuations of magnetic ions. There 
the role of heavy particles was played by magnons; com- 
plete spin polarization of conduction electrons prevented 
magnon emission or absorption processes, and only the 
processes of electron-magnon scattering being allowed. 
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